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Abstract 

Taking into account the occurrence of a zero of the surface diffusion current and 
the requirement of the Ehrlich-Schwoebel effect, Siegert et al [23j formulate a model 
of Langcvin type that describes the growth of pyramidlike structures on a surface 
under conditions of molecular beam epitaxy, and that the slope of these pyramids 
is selected by the crystalline symmetries of the growing film. In this article, the 
existence and uniqueness of weak solution to an initial boundary value problem for 
this model is proved, in the case that the noise is neglected. The regularity of the 
weak solution to models, with/without slope selection, is also investigated. 
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1 Introduction 

Many processes occur mainly at surfaces of materials, such processes include crystal 
growth, catalytic reactions, production of nano-structures. Thus surfaces are of great 
technological and fundamental interest. There has been increasing interest in the un- 
derstanding of the kinetics of surface growth processes, see e.g. |14l [2]. The growth of a 
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crystalline film from a molecular or atomic beam is commonly referred to as molecular 
beam epitaxy (MBE), which is among the most fined methods for the growth of thin 
solid films and is of great importance for applied studies [9]. MBE takes place in high 
vacuum or ultra high vacuum, for instance, some 10~* Pa. The most important aspect 
of MBE is the slow deposition rate that is typically less than lOOOnm per hour, so MBE 
allows the films to grow epitaxially. The slow deposition rates require proportionally 
better vacuum to achieve the same impurity levels as other deposition techniques. In 
turn, it is possible by using this technique to grow high-quality crystalline materials and 
form structures with very high precision in the vertical direction. 

The mathematical theory has been developed since Langevin who proposed the equa- 
tion named after him. The model by Kardar, Parisi abd Zhang (KPZ) |10j describes 
well growth process such as the Eden process [5j, ballistic deposition [26], and growth 
of various restricted solid-on-solid models [19]. Thus the KPZ model has been widely 
accepted as a model for the growth of crystals, and also has been extended to various 
cases. Let us mention especially the following work: In [9] the authors consider the sys- 
tem that has potential barriers near step edges that suppress the diffusion of adatoms 
to a lower terrace. This effect is now commonly called Ehrlich-Schwoebel effect. They 
take into account the step-flow regime and instability, and proposed a continuum equa- 
tion to model the growth in MBE, which is valid only at the early time as long as the 
slopes are much smaller than 1. However this is too restrictive to describe the unstable 
three-dimensional growth of real materials in the later time regime. In order to take into 
account the occurrence of a zero in the surface diffusion current and the requirement of 
Ehrlich-Schwoebel effect, Seigert et al |23] thus introduce a current for a structure with 
a cubic symmetry so that the current is changed to the one that has a zero differing 
from so that the model can be also applied to the regime when the slope is much 
greater than 1. In this article, we shall study this model, also the model without slope 
selection will also be investigated. Due to the forming of steps, pyramidlike surfaces, 
etc. during the crystal growth, it is more natural to assume that the initial data is in 
L'^{n) or H'^{Q) than in H^{n). 

To formulate the model, we need some notations. Let x = (xi,X2) G ri C be 
material point, is an open bounded set with smooth boundary dQ. Let t be the time 
variable. Qt = (0, t) x 0,. h = h{t, x) is the height that is measured in a co-moving 
frame of reference and describes the local position of the moving surface. V^^/i is the 
gradient of h, and A^,. = d"^ /dx\ + d"^ ldx\ is the Laplacian. For simplicity of notations, 
we shall use the following notations 



which is satisfied in Qt; where T is a given positive number. And the boundary and 
initial conditions are 



Then the equation turns out to be 




+ vl^lh + div^ {J{Vxh)) = 0. 



(1.1) 



d_ 
dn 



h = 0, {uV^A^h + J{Vxh)) ■n = Q 



on [0, T] X an, 



(1.2) 
(1.3) 
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Here n = (ni,n2) is the unit outward normal vector to the boundary dQ. We have 
introduced the surface diffusion current J = J{Vxh) = J{p, q) 



J = (ji, i2), (1.4) 



Jl = 

h = 

where / is defined by 



«((P + q)f {{P + q?) + {P- q)f {{p - qf) ) , (1-5) 
a{{p + q)f{{p + qf)-{p-q)f{{p-qf)), (1-6) 



and a is a constant of surface diffusion, /3 = (id)'^ where id is the diffusion length. 

This completes the formulation of an initial-boundary value problem. It is worth a 
remark on the nonlinearity / since there are several varieties of / which lead to different 
models relating to ours. 



Remark 1. Define 



Then the corresponding model is proposed by Johnson et al. [9j . However this model does 
have a slope selection mechanism, and is correct only for early time as long as the slopes 
are much smaller than 1. It is too restrictive to describe the unstable three-dimensional 



growth of real materials in the later time regime. Therefore, (1.7) is introduced by Siegert 
et al in to interpolate the two regimes. The form of the surface current J is the 
minimal model in the sense that the nonlinearity must be chosen such that it describes 
the instability and leads to slope selection. The flux still has the correct physical behavior: 
\J\ ~ Y^p2 _|_ q2 Jqj. p2 _|_ ^2 ^ _L and \ J\ ~ 1/ \/p^ + for ^ <C + (7^ ^ 1 . This type 



of fluxes gives rise a completely different behavior than the one defined by (1.8), despite 
many similarities, as shown in [23j . The exact form of f does not play a role since the 
slope selection mechanism and the growth exponents do not depend on such details. 
In |241 120| . the current of the form 

J = e(l-|eP) (1.9) 

is used, however, it has stable zeros for all slopes with |^| = 1 regardless of the direction 
of S,. Thus such an azimuthal symmetry is unrealistic for crystalline films. Therefore 
here fi (i = 1,2) are functions chosen such that fi{p'^,q'^) = f2{q^-,p'^)- The simplest 
form that describes growth on such substrates is a current with components 

jl = p{l-p^-bq^), (1.10) 
j2 = q{l-q^-bp^), (1.11) 

which leads to a buildup of pyramids with selected slopes {po,qo) = (±1, ibl)/\/l + b for 
— 1 < 6 < 1. This diffusion current is suitable for substrates with a quadratic symmetry. 
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Finally we point out that after a coordinate transformation X = Ax, where 



we find J in (1.7) can he reduced, without loss of generality, to a simpler form 

h = apfip^), (1.12) 
J2 = aqfiq^). (1.13) 

It is interesting to compare the important difference between the Cahn-Hilliard equa- 
tion modehng phase-ordering and the model considered here. MathematicaUy, equation 



(1.1), with nonhnearity (1.8), (1.9) or (1.13), differs from the Cahn-Hilhard equation 



due to the flux term: J in this paper depends on the gradient of the unknown, while 
the Cahn-Hilliard equation on the unknown only. Numerical experiments also show 
the important differences between this model and the Cahn-Hilliard one. Many papers, 
e.g. [23l [22l [21] have been carried out the study of the differences between the problem 
studied here and phase-ordering dynamics described by the Cahn-Hilliard equation [3j. 
These differences become apparent when the domain configurations are plotted as in 
Figure 1. A domain in this context is an area of constant slope corresponding to one of 
the four values. The analogous case in phase-ordering dynamics is described by a four- 
state clock model, see, e.g. |2H[8]. However, in that case we shall find that domain walls 
do not have any particular orientation, whereas here domain walls are intersections of 
planes of constant slopes and therefore form straight lines. Furthermore, there are two 
types of domain walls: Domain walls at which only one component of the slope changes 
are aligned along the x and y axes. These are the edges of the pyramids; Domain walls 
at which both components of the slope change run at 45° with respect to the principal 
axes. These latter domain walls form the so-called roof tops. 




Figure 1. Left: the configuration for the crystal growth model, from [22j; Right: the 
domain walls of phase-ordering governed by the Cahn-Hilliard equation, from e.g. [8]. 

Now we derive the model briefly. Define the free energy by 

E[h] = ^ (|(A,/i)2 - I (F {{hx, + hx,f) + F {{hx, - K,f))) dx. 

Here F = fiy)dy is the primitive of / such that F(x) grows at — log(|a;|) as |x| goes 
to oo. Suppose now that h = h{t,x) is a solution to ( |1.1| ) - (1.3). Formal computations 
yield 



dt 



E[h] 



(^-af {{K^ + K^f) {hx-, + K.,){hx^ + hx2)?jdx 



+ j i^-af - hx2f) {hxi - /ix2)(^xi - hx2)t + i^A^hAxht^dx. 

Combining the terms containing {hxj^)t (or {hx2)t) together we can rewrite the right 
hand side of the above equality as 

"(/((^xi + hx2f){hx^ + hx2) + / ((^xi - hx2f) {hx^ - hx2)^ {hx^)t dx 

"(/ ((^2^1 + hx2)'^) (hxi + hx2) - f {{hxi - hx2f) {hxi - hxifj ihx2)t dx 

+ / vAxhAxhtdx 

J ■Vxht + vAxhAxht)dx. (1.14) 



Using integration by parts and equation (1.1), we infer from (1.14) that 
d 



dt 



m 



/ {—J — i^^xAxh)-'Vxhtdx 
Jn 

/ divx {J + i^^xAxh) ht dx 
Jn 

- / {divxiJ + i^yxAxh))^ dx 



< 0. 



(1.15) 



This implies the second law of thermodynamics is valid. The equation can be written 
in a gradient form with the total surface current defined by Ji := Jeq + JiVxh) where 
Jeq := v^xAxh is called the equilibrium surface current, v is proportional to the surface 
stiffness. 

Statement of the main result. Before the statement of our main results, we define 
weak solutions to problem consisting of ( 1.1 )-( 1.3). We use the notations: {f,g)Qj, and 
{f,g)n are, respectively, the inner product of / and g over Qt and Q. {f,g)Qj, denotes 
the dual product of f,g with / G L'^{0,T; X'), g £ L'^{0,T;X) and X is a Banach 
space and X' is its dual. H"^(Q) are the standard Hilbert spaces of order m. Define 
H%{n) := {/ G H^{n) I ^ = on dn} and its dual space is denoted by H'^iQ). 
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Definition 1.1 Let ho € L^(il). A function h = h{t,x) with 



h G L^{0,T;L'^{n))nL'^{0,T;H^{n)), ht G L\0,T; H]^^{n)) 



(1.16) 



is a weak solution of problem ([Op - ( fj.3| j, if 

{ht, ip)Q^ + i/(A^/i, A^(/?)q^ + (J, V^99)q^ = (1.17) 

holds for all G L^(0, T; //^(fi)), and limt^Q+ (/i(t), -0)^ = (/io,V')n /or ^ L'^{^)- 

Now we are in a position to state the main results of this article. 

Theorem 1.2 (Existence) Suppose that the boundary of is smooth, and ho G 
L^(r2). Then there exists a unique weak solution h of problem (1.1) - (1.3) in the 
sense of Definition 1.1, and the total mass j^h{t,x)dx is conserved, i.e. J^h{t, x)dx = 
ho{x)dx. 

Moreover, if h^ G H^{il,), the weak solution of problem (1.1) - (1.3) has, which in 
addition to (1.16), the following regularities 

he L°^{0,T;H\n))nL\0,T;H^{n)), ht e L'^{0,T; H-\n)). (1.18) 

Theorem 1.3 (Regularity) Suppose that ho G H'^"''{Q) with m G N. Then the weak 
solution h of problem ( |1.1[ ) - (1.3) satisfies 



h G L°^(0, T; H^'^iQ)) n L^(0, T; F^™+^(f))), ht G L^(0, T; //^'""^(rj)), (1.19) 



r2m+2 



2m-2/ 



and 



D^/i G L°°(0, T; L2(r?)) n L'^{0, T; //^(r?)), if m = 2l, / G N; 

G L°°(0, T; H^{n)), D[h G ^^(o, T; L'^{n)) if m = 21 -I, / G N. 

Consequently, if G C°°(r2), i/ien i/ie solution h is smooth on Qt. 



(1.20) 



Now let us recall some references related closely to our problem. In |llj an initial 
boundary value problem of this epitaxial model with cubic nonlinearities is studied in 
which the initial data is chosen in H'^{Q), and the i?^(ri)-norm of the solution follows 
directly from the Clausius-Duhem inequality, the second laws of thermodynamics, but 
this technique does not work for our case since we assume the initial data is only in 
H^[Q). Li and Liu study the initial boundary value problem for the MBE model with or 
without slope selection in [T7j and the boundary conditions are chosen periodic. In both 
articles, they construct approximate solutions by using the Galerkin method, while we 
use a linearized problem, together with the convolution technique, to obtain a sequence 
of smooth approximate solutions, then establish a priori estimates for this sequence. 
Kohn and Otto |12j investigate the coarsening Rate for the Cahn-Hilliard equation, and 
Kohn and Yan |13j studies the coarsening rate for an epitaxial growth model with a 
cubic nonlinearity. Watson and Norris [27] study the coarsening dynamics of multiscale 
solutions to a dissipative singularly perturbed partial differential equation with a trigo- 
nally symmetric potential which models the evolution of a thermodynamically unstable 
crystalline surface. 
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We give now a remark on the choice of the initial data in this article. 
Remark 2. The assumption that initial data Hq is in is more natural than the one 
that ho is in . The reason is that h is piecewise affine in the case that the surfaces are 
high-symmetric, such as pyramidlike ones. One evidence can be also seen from a typical 
Scanning Tunneling Microscope (STM) picture, see Figure 2, which shows clearly that 
the surface are not smooth. Correspondingly, a good mathematical model should consider 
this feature. 



If the initial data is in , the existence of weak solutions to problem (1.1) - (1.3) 



in which the nonlinearity is cubic, like (1.9), also (1.10) - ( jl.lip , is still open, and may 



be interesting. To solve such a problem, I surmise we need to invent an inequality of 
the Brezis-Gallouet type [3j. 




Figure 2. Typical STM image of mounds appearing on the surface. From Krug, Politi 
and Michley [15] . 



The organization of the remaining parts of this article is as follows. The existence and 
uniqueness of weak solutions is studied in Section 2 by constructing smooth approximate 
solutions and using a priori estimates. Section 3 consists of two parts, each of which 
is concerned, respectively, with the regularities of weak solutions to models with and 
without slope selection. 



2 Existence of weak solutions 

2.1 Existence for the approximate problem 

To prove the existence Theorem |1.2[ we construct smooth approximate solutions in 
a similar way as done in \T], provided the initial data is smooth. Then we establish 
uniform a priori estimates of these solutions by which we conclude compactness. Before 
formulating an approximate problem, we introduce the modifier x = xit^^) such that 
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X G C^{Qt) satisfies x{'^-,x)dtdx = 1. For e > 0, we set 

Xe{t,x) 

and for any function / G L°°{Qt) we define 



1 t X 

£ \e e 



f{t,x) = ixe* f){t,x) = / Xe{t- s,x-y)f{s,y)dsdy. 



(2.1) 



We choose a smooth sequence such that 

11^0 ~ ^ollx - 

as e ^ 0. Here X = L'^{n) or X = H^n). 

Then the smoothed problem turns out to be 



dh 



+ i/A2/i + div^(J(V^/i)) = 0, 



with the following boundary and initial conditions 
d 



dn 



h 



0, ( vV^^^h + J{V^h)j - 71 = 0, 
/i(0,x) = /ig(x), 



on [0,T] X d^, 
x gQ. 



(2.2) 

(2.3) 
(2.4) 



Note that equation (2.2) is a linear fourth order parabolic equation with a smooth 



known term. By the existence theorem for higher order parabolic equations in the book 
by Ladysenskaya et al. jH] or the book by Eidelman [6], we have 

Theorem 2.1 (Existe nce of smooth approximate solutions) Suppose that the as- 
sumptions of Theorem 1.2 are satisfied, and e is a given positive constant. Let h € 
L'^{G,T-H^{n)). 



Then for any T > 0, there exists a unique smooth solution h of problem (2.2) - (2.4), 
which satisfies that the total mass h{t,x)dx is conserved. 

The solution h constructed in Theorem 2.1 depends on the small parameter e. In 



order to prove the existence of weak solutions to the original problem (1.1) - (2.4), we 



need to establish some a priori estimates which are independent of h and e and thus 
guarantee the passage to limit of as e ^ 0. 

2.2 A priori estimates 



Assume that there exists a classical solution /i^ to problem (2.2) ~ (2.4) with smooth 
initial data /iq and h G L'^{0, T; H'^{Q,)) satisfying 



lL2(o,T;_f/2(r2)) ^ C. In what follows, 
C denotes a constant which is independent of e and h. The L'^(il)-norm of / is denoted 
by 11/11- We shall derive a priori estimates for this solution. To begin with, we first state 
the following lemma on the nonlinearity /. 
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Lemma 2.1 There hold, for all y > 0, that 

\yf{y)\ < c, (2.5) 
< C,meN. (2.6) 

Here /('"^ is the m-th order derivative of f . 

Proof. One needs to investigate the behavior of the function 

-yf{y) ^ ^ 



(l-y)2+/3y 

for all y > 0. It is easy to see that —yf{y) — )• 1 as y — )• oo, hence | — yf{y)\ < C 
for y > M, where M is a suitably large constant. We shall prove that —yf{y) is also 
bounded on the interval [0,M]. To this end, we need to prove that the denominator 
g{y) = (1 — y)^ + (3y is greater than a positive constant. Note that g{y) is nonnegative 
for y > 0. Thus is the only possible minimum of g{y). However g{y) = implies that 
1 — y = and y = which cannot be satisfied simultaneously. Therefore, the minimum 
of g must be positive, and since M is finite, we infer from the continuity of g{y) that 

min fl(y) = Ci > 0. (2.7) 

0<y<M 

We then obtain, for ah y G [0,M], that 

\yf{y)\ < < C. 

Next we consider the behavior of derivatives of /. Let m < 1 be an integer. Rewrite 
where g{Y) = and Y = {1 — y)^ + /3y. Invoking the product rule 



m 

/r^=(/i-/2)(™)=E^-/^'v^"'^ 

k=0 



where denotes the number of /c-combinations of an m-element set, we have 

/(™)(y) = C^(l-y)(y(y(y)))(™) + Ci(l-y)'(y(y(y)))(™-i) 

= {l-y)igiY{y)))("^^-migiYiy)))("^-'\ (2.8) 

Making use of the Faa di Bruno formula, i.e. 

and recalling (^(y))^™) = (y-i)('^) = (-l)™m! y-^-i, Y'{y) = -2(1 - y) + /3, 
Y"{y) = 2 and Y^^^y) = for any j > 2, we can reduce (|2.9l) to 

E (-i)'-'=^^r-".--"te)n(^)''.(2.io) 

lh+2l2=m j=l\ / 
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For y E [0,M], (y) is smooth by (2.7). One thus needs only to investigate the 
behavior for large y G [M, oo), and it is enough to calculate the highest exponent. From 
(2.10) it follows that the highest exponent is less than or equal to 

-2(^1 + /2 + 1) + l/i + 0/2 = -h - 2/2 - 2 = -m - 2. 



Therefore, invoking (2.8), we assert that there exists a constant 7 such that 



(2.11) 



Hence, this implies (2.6). Thus the proof of this lemma is complete. 

From now on we are going to derive a priori estimates. The first is 
Lemma 2.2 There hold for all t G [0, T] 



h{t,x)dx = / hQ{x)dx, 
Jn 

+ / \\AMT)fdT < C. 

Jo 



Proof. Integrating (|2.2|) with respect to x yields 

h{t, x)dx = 0, 



d^ 
dt 



(2.12) 
(2.13) 

(2.14) 



which implies (2.12). 

Multiplying (2.2) by h and integrating the resulting equation with respect to x, we 
obtain 



+ iy\\A^hf+ / J{Vcch) -V^hdx = 0. 



ld_ 

Applying Lemma 2.1 we obtain 

J{Vxh) ■ Vxhdx < C\\Vxh\\. 



(2.15) 



(2.16) 



Using the Poincare inequality of the form ||Va;/i|| < C||D^./i|| + C{f^ hdx), applying the 
elliptic estimates < CHA^^/iH and from ( |2.12 ) it then follows that 

\\V^h\\ < C\\A^h\\ +C. 



Therefore, (2.15) becomes 

1 d 
2dt 



\hf + u\\A^hf < -\\A^hf + C. 



(2.17) 



From this, estimate (2.13) follows. And the proof of this lemma is complete. 
Lemma 2.3 There holds 



\\HtWmin)+ / Wy.AMrWdr < C. 



(2.18) 
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Proof. Integrating, with respect to x over 17, equation (2.2) multiplied by — A^/i and 
using integration by parts yield 



\ d f ' — ' — " 



Making use of Lemma 2.1 again, we can easily prove that 

J(V,/i) • V,A.,hdx < C + ^||V,A,/if. 



n 



(2.19) 



(2.20) 



Thus integrating (|2.19|) with respect to t one gets 

r-t 



^ Jo ^ JO 



(2.211 



From which estimate (2.18) follows. The proof of this lemma is complete. 



2.3 Existence of solutions to the phase field model 

In this section we shall make use of the a priori estimates established in the previous 
subsection to study the convergence of the solutions h"^ of the approximate problem for 



m — >■ oo, thereby proving Theorem 1.2 In our investigation we need the following well 
known results, see, for instance. Lions [18j, Evans [7]: 

We shall make use of the following lemma which is of Aubin-Lions type. 

Lemma 2.4 Let Bq, B, Bi be Banach spaces which satisfy that Bq, Bi are reflexive 
and that 

BqCCB C Bi. 
Here, by CC we denote the compact imbedding. Define 



W 



f\feL^{0,T;Bo), |Gr(0,T;B; 



with T being a given positive number and 1 < r < oo. 
Then the embedding ofW in C([0,T];i?) is compact. 

To deal with the nonlinear terms, we also need 

Lemma 2.5 Let T be an open set in . Suppose functions fn, f are in L'^(T) for any 
given 1 < q < oo, which satisfy 

ll/n||L'j(r) ^ C) fn^f almost everywhere in T. 
Then fn converges to f weakly in L^(r). 

We now turn to prove the existence of weak solutions for the initial /iq that is assumed 



in L^(J7). The existence for the initial data in H^{Q,) is easy by recalling Lemma 2.3 



Using problem (2.2) - (2.4), we can construct smooth approximate solutions as follows: 



Let h^ be a given function. By solving problem (2.2) - (2.4) with h = h^, one gets h^. 
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Suppose we have obtained /i™ for some m G N, then set e = ^ and h = h"^. Thus, 
we can define /i™'"*"^ successively. Therefore we has a sequence of approximate solutions 
/i*". Since equation (1.1) is nonhnear, we need some results about strong and pointwise 
convergence. 

Define f = and r = 2. Set 

Bo = H^in), B = H^in), Bi = H-^{n), 

it follows from Lemma 12.21 that 

Il2(o,T;Bo) - ^' \\ht\\L'2{0,T;Bi) < C. 



Applying Lemma 2.4, we then conclude that {h"^} is compact in C{[0,T]; B), namely 
C{[0,T];H^{n)), so {/i^}, for i = 1,2, is compact in C{[0,T]; L'^{^})). Therefore there 
exists a subsequence, still denote it by {/i™}, such that 







(2.22) 



as m — )■ oo. Moreover, we can select furthermore a subsequence such that /i^ converges 
to hx^ almost everywhere. Setting = ^, 5"' = Vxh"^. By the properties of convolution 
and (2.22), we have 

\\xn*s--su2(Q,j < \\x.*{s''-s)\\mQ,j + \\{s-x.*s)\\mQ,^) 

< \\(S-x.*S)h2^Q^j + \\S'^-Sh2^Q^j^0, (2.23) 

for K — )• 0, whence we can select a subsequence, still denote by X/t * S'^ converges to S 
almost everywhere. Consequently, we assert that 

Jiy xh"^) converges to J{Vxh) 

almost everywhere as m — oo. Remembering that \\JC^xh"^)\\L^{QT) < C, using 
Lemma 12.51 we assert that 



J(V^/i^ 



J{Vxh) 



in LP'{Qt) as m — ^ oo. 

For the linear terms, by weak compactness, one can easily get 

{hf,^) {ht,ip), (Axh"", Ax^) (Axh, Axip) 

as m oo, for ah ip G L'^{0,T; Hf^{0,)). 



(2.24) 



(2.25) 



Taking the inner product of (2.2) and if we arrive at 

= { — ,^) + u{Axh"',Ax^)-{JiVxh"'),Vx^) 



dt 
dh 

( — , 99) + u{Axh, Axip) - {J(yxh),Vx^). 



(2.26) 



Thus (1.17) is proved. From ( 2.12| ), (2.22) and the choice of the smooth initial data h] 
(let e = — ) we have 



h"^{t,x)dx = / h^{x)dx 



ho{x)dx, 



(2.27) 
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and the left hand side converges to J^h{t, x)dx, thus the mass is conserved for weak 
solution. And the existence of weak solutions is complete. 

Next we are going to study the 
Stability and Uniqueness. Let hi, /i2 be two weak solutions corresponding to initial 
data /iq and /ig, respectively. Define u = hi — h2. We write 

^xhi = {pi,qi), {i = l,2), VxU={p,q), Ji = J{Vxhi). 



Then by the estimates in Lemma 2.1 we have | Ji — J2I ^ C\^xhi — Vxh2\ = C'|Va;ti| 
hence 

= - — ||n||^ + z^llArnll^ + / ( Jt — Jo) • V rudx 



2 dt 



\uf + i^WAxuf + / (Ji-J2)-V, 
Jn 



> ^^\\uf + i^WAxuf - C\\Vxuf- (2.28) 
Making use of the Nirenberg inequality of the following form 

||V:^n|| < CIIA^iillsllnlls +C"||n||, (2.29) 



and the Young inequality, from (2.28) one obtains 

~\\uf + i^WAxuf < C(||A^u||||n|| + ||nf ) 

< '^WAxuf + C\\uf . (2.30) 

Now using the Gronwall inequality we get 

\Ht)f < ||u(0)||2e^*. (2.31) 

Here u{0) = /ij — /iq. Thus the solution depends continuously on the initial data. 

Consequently, if Hq = hQ, that is ||n(t)|p = which implies ||M(t)|p = 0, so the weak 
solution is unique. Therefore, the proof of Theorem 1.3 is complete. 



3 Regularity of weak solutions 

We shall investigate the regularity of weak solutions to both models with and without 
slope selection, while for the latter model we can only carry out such study in one space 
dimension. 



3.1 The model with slope selection 



Suppose now that ho € H {Q) with m E N and there exists a unique solution h to 
problem ( 1.1 ) - ( 1.3 ). In this section we shall investigate the regularities of this solution. 
We first consider the case that m = 1. 

Lemma 3.1 There hold for ho £ //^(Jl) that 



+ \\h 
\\h\ 



*IIl2(0,T;L2(q) 
2 

L2(o,r;H4(n)) 



< 
< 



c, 

C. 



(3.1) 
(3.2) 
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Proof. Multiplying equation (1.1) by ht and integrating the resulting equation with 
respect to x yield 







\htf + ~\\A^hf + I J-V^htdx 
\htf + ^^E[h]{t). 



Thus one has 

/ Whfdt < c, 

Jo 

E[h]{t) < E[h]{t). 
Recalling the definitions of E[h] and / we arrive at (|3.1 



(3.3) 



From equation (1.1), making use the estimates in Lemma 2.1 and (3.1) we obtain 

[ \\Alh{T)\\'^dT < C [ \\ht{T)fdT + C [ WD^hfdT, 
Jo Jo Jo 

< C. 



Thus (3.2) is proved. And the proof of this lemma is complete. 

To get the a priori estimates for higher order derivatives, we differentiate the equation 
with respect to t to get 



dt 



+ vAlht + div^ (( J(V^/i))t) = 0. 



(3.4) 



Such computations are formal. However by using the technique of difference quotient 
one can justify easily. In a similar way for deriving the estimates for /i, we arrive at 

Lemma 3.2 Suppose that /iq G H^{VL), i.e. m = 2. There holds for any t G [0, T] that 



\htr + 



H4(o) + ll^i|lL2(0,T;H2{Q) - ^■ 



(3.5) 



Proof. Multiplying equation (3.4) by ht and integrating the resulting equation with 
respect to x yield 







^j^Whf + i^WA^htf + I JfVMx 

1 d 
2~dr 



n 



> --\\htf + u\\A,htf - C\\VM\\ 



(3.6) 



With the help of the Nirenberg inequality (2.29) one obtains from (3.6) and the 
Young inequality that 



1 d 
2dt 



\htf + u\\A^htf < C\\A^ht\\\\ht\\+C\\ht 
< "^WA^htf + CWhf. 



(3.7) 
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From which, by the Gronwah inequahty, it fohows that 

rT 



\h. 



AM-^dt < C. 



(3.8) 







Furthermore, one can get from equation (1.1) and Lemma 2.1 that 

|2 



\Alh\ 



< C\\ht\ 

< C, 



+ C\\Aih\\^ 



(3.9) 



which implies that h £ L°°(0, T; H'^^Q,)). The proof of this lemma is thus complete. 
For the higher order derivatives with both x and t, we have 

Lemma 3.3 Let ho £ H'^"^{n). There holds for any t G [0,T] that 

||/l(t)||^2m(Q) + ll^llL2(0,T;_H"2™+2(n) + ll^*llL2(0,T;H2™-2(n) - ^' 



(3.10) 



and 



\Dih\\ 



L°°(0,T;L2(q)) 



+ \\DiHL2(0,T;H'^(n) 



< C, if m = 21, (3.11) 

< C, if m = 2l-l. (3.12) 

Proof. We employ the mathematical induction. From Lemma 3.1, it is easy to see that 
and (3.12) are true for m = 1 which implies I = 1 too. By Lemma 3.2, estimate 



3.10 



3.11 



3.12 



holds when m = 2. Assume that (3.10) is true for any A: < m G N and (3.11 ) and 
are true respectively for even and odd m. Next we shall prove they are true for 



k<m + l when ho € 



m+l 



(resp. 



times 



For the case that m + 1 is even (resp. odd), differentiating 

m-f-l m 

equation (1.1) with respect to t, letting v = D^'^ h (resp. v = h), repeating the 
argument of Lemma 3.2 (resp. Lemma 3.1) for this function v, and using the estimates 
in Lemma 2.1, we then conclude (3.10) holds for m + l, moreover, (3.11) (resp. (3.12)) 
is true for / = (resp. / = y). Thus the proof of this lemma is complete. 

3.2 The one-dimensional cubic model without slope selection 

In this subsection we shall study the one-dimensional problem with a cubic current, 
i.e. ( |1.9[ ), or ( |1.10[ ) - ( |1.11[ ), which is studied in [Ml EOj [27]. However the original two 
dimensional problem is still open. The problem is 



ht + i^h^ 



XX 

hx — 0, hxxx 
h\t=o 



in Qt, 

on [0, T] X dn, 

ho- 



(3.13) 
(3.14) 
(3.15) 



Here ^2 = (a, 6) C M, a, 6 G M and J = ahx{l - hi). 

Assume that ho G L^(il). Multiplying ( [3.13 ) by h and integrating it with respect to 
X give 

1 d 

2dt' 







1 II, ii2 II, ii2 

2dt " 



J{hx)hxdx 
\hx\'^dx — \\hx 



(3.16) 
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By the Young inequality < ia^ + i, we infer from (3.16) that 



2dt 
which gives 

Lemma 3.4 We have 



+ J^\\hxx\\'^+ [ \hx\'^dx<}- [ \hx\'^dx+'^, (3.17) 



+ 



2 



+ / \hXdx \dt<C. 



(3.18) 



Based on this lemma we can define weak solutions and prove the existence and 
uniqueness in a similar way as in Section 2 for /iq G L^($7). Suppose that /iq G i?^($7), 
is the weak solution regular, i.e. h G L^(0, T; H^(Q,))7 To answer this question, we need 
more estimates. From equation (3.13), we obtain 

1 d 







Therefore, 



2dt 
1 d_ 
2dt 



+ 



\h ||2 _|_ jAlh ||2 
\"'x\\ 1^ \ \' ''XxxW 



n 



J ( )xhxx dx , 



n 



+ / 3|/ixl l/ixxl ] dt<C. 



(3.19) 
(3.20) 



So the weak solution is more regular if ho G H^{U). 

However for two-dimensional problem, we can not get (3.19) so that we obtain the 
good term f^Slhxl'^lhxxl'^dx in (3.20). 
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